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Abstract 

We investigate the validity of the generalized second law (GSL) of gravitational thermo- 
dynamics in a non-flat FRW universe containing the interacting generalized Chaplygin gas 
with the baryonic matter. The dynamical apparent horizon is assumed to be the boundary 
of the universe. We show that for the interacting generalized Chaplygin gas as a unified 
candidate for dark matter (DM) and dark energy (DE), the equation of state parameter can 
cross the phantom divide. We also present that for the selected model under thermal equilib- 
rium with the Hawking radiation, the GSL is always satisfied throughout the history of the 
universe for any spatial curvature, independently of the equation of state of the interacting 
generalized Chaplygin gas model. 
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1 Introduction 



One of interesting DE models is the Chaplygin gas (Kamenshchik et al. 2000, 2001; Bilic 
et al. 2002) and its generalization (Bento et al. 2002, 2004) which has been widely studied 
for interpreting the accelerating universe. It is remarkable that the generalized Chaplygin gas 
(GCG) equation of state has a well defined connection with string theory and can be obtained 
from the light cone parameterization of the Nambu-Goto action, associated with a D-brane (Bilic 
et al. 2007). The GCG is the only gas known to admit a supersymmetric generalization (Zhang 
et al. 2006). In this model, a single self-interacting scalar field is responsible for both DE and 
DM, giving also the observed accelerated expansion (Garcia-Compean et al. 2008). The striking 
feature of the GCG is that it allows for a unification of DE and DM. This point can be easily 
seen from the fact that the GCG behaves as a dust-like matter at early times and behaves like 
a cosmological constant at late stage. This dual role is at the heart of the surprising properties 
of the GCG model (Wu and Yu 2007; Lu et al. 2009). Moreover, the GCG model has been 
successfully confronted with various phenomenological tests involving SNe la data, CMB peak 
locations, gravitational lensing and other observational data (Alcaniz et al. 2003; Bento et al. 
2003a,b; Bertolami et al. 2004; Bento et al. 2005; Wu and Yu 2007; Lu et al. 2008; Lu et al. 
2009). 

Here our aim is to investigate the GSL of thermodynamics for the GCG model as a candidate 
for the unified DM-DE which is in interaction with the baryonic matter in the non-flat universe 
enclosed by the apparent horizon. 



2 Interacting GCG and BM 

The GCG model is based on the equation of state 

^Ch = ~ (1) 

where A and a are the GCG constant parameters (Bento et al. 2002, 2004). The case a = 1 
corresponds to the standard Chaplygin gas model (Kamenshchik et al. 2000, 2001; Bilic et al. 
2002). 

Using Eq. (1), the continuity equation can be integrated to give 



Pch = PChO 



A s + 



[I- A 



i 



s 



a 3(l+a) 



l+a 

(2) 



where A s = -i^, /?cho is the present energy density of the GCG, and a is the cosmic scale 

^ChO 

factor. Note that the GCG model smoothly interpolates between a non-relativistic matter 
phase (pch a~ 3 ) in the past and a negative-pressure DE regime (pch = ~Pch) at late times. 
This interesting feature leads to the GCG model being proposed as a candidate for the unified 
DM-DE (UDME) scenario (Wu and Yu 2007; Lu et al. 2009). 

In the framework of the Friedmann-Robertson- Walker (FRW) metric 

ds 2 = _ dt 2 + a 2 {t) + r 2 dQ 2\ ? (3) 



the first Friedmann equation has the following form 

8vr 

f / 2 ~3 



tf 2 + 3 = — (Pch + Pb), (4) 
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where we take G = 1 and k = 0, 1, —1 represent a flat, closed and open FRW universe, respec- 
tively. Also pet and pb are the energy density of GCG and BM, respectively. In terms of the 
the dimensionless energy densities as 

^ b= 3i^' ^ Ch = ^i^' nk = ^lP> (5) 
the first Friedmann equation yields 

n h + n Cil = i + fifc. (6) 

The energy equations for the GCG and BM with uj^ = are 

PCh + 3#(1 + w C h)pch = -Q, (7) 

p b + 3Fp b = Q, (8) 

where following Kim et al. (2006), we choose Q = 3b 2 H(pch + pb)- The interaction between the 
GCG and BM can be detected by local gravity measurements (Hagiwara et al. 2002; Peebles 
and Ratra 2003). Following Guendelman and Kaganovich (2008), the interaction between the 
BM and the DE will cause the appearance of the fifth force. However, according to the results 
of the fifth force experiments, a coupling of the DE (modeled by a light scalar field) to visible 
(baryonic) matter is strongly suppressed. According to Guendelman and Kaganovich (2008), a 
new aspect introduced by modern cosmology to this problem is the question of why the coupling 
of the light scalar (DE) to visible matter is strongly suppressed while similar coupling to DM 
is energetic. Discovery of DE and cosmic coincidence interpreted as evidence of the existence 
of an unsuppressed DE to DM coupling, turns the fifth force problem into an actual and even 
burning fundamental puzzle. 

Taking time derivative of Eq. (2) yields 



o(i a \ u Pew 
Pch = -3(1 - A s )H-y——j 



a 3(l+a) 



l+a 

(9) 



Substituting Eq. (9) in (7) gives the equation of state (EoS) parameter of the interacting GCG 
model as 



Observational constraints on the GCG model as the UDME from the joint analysis of the latest 
astronomical data give the best-fit values of the GCG model parameters. In this respect, Wu and 
Yu (2007) obtained (A s = 0.76, a = -0.005) for non-flat universe with (Q k = 0.04, ft Ch = 0.995). 
Also Lu et al. (2009) derived (A s = 0.73, a = -0.09) for flat universe with (fi C h = 0.956). 
Equation (10) shows that in the absence of interaction between GCG and BM, b 2 = 0, and for 
the present time, a = 1, uch = —A s > — 1 and cannot cross the phantom divide. However, in 
the presence of interaction, b 2 ^ 0, taking (A s = 0.76, a = -0.005, n k = 0.04, Q Ch = 0.995) 
given by Wu and Yu (2007), (A s = 0.73, a = -0.09, fl k = 0.0, ft C h = 0.956) given by Lu et al. 
(2009) and a = 1 for the present time, Eq. (10) clears that the phantom EoS, i.e. ojQh < — 1> 
can be achieved when b 2 > 0.23 and b 2 > 0.26 for non-flat and flat universe, respectively. 
The deceleration parameter is given by 

H 



(id 
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Taking time derivative of both sides of Eq. (4), and using Eqs. (5), (6), (7) and (8), we get 



Ui + n k + 3n Ch Lo Ch ). (12) 



Substituting Eq. (10) in (12) yields 

-3^,a 3 ( 1+Q ) 



q= 2 



1 - A 8 + A.a^+a) 



n ch + (i-3b 2 )(i + n k ) . (13) 



Now taking (A s = 0.76, a = -0.005, n k = 0.04, C h = 0.995) given by Wu and Yu (2007), 
(A s = 0.73, a = -0.09, tt k = 0.0, Q C h = 0.956) given by Lu et al. (2009) and a = 1 for the 
present time we get q = —0.61 — 1.566 2 and q = —0.55 — 1.5b 2 for non-fiat and flat universe, 
respectively. These results show that the deceleration parameter is always negative even in the 
absence of interaction between GCG and BM. Therefore the GCG model in the present time 
can drive the universe to accelerated expansion. 



3 GSL of thermodynamics 

Here, we study the validity of the GSL in which the entropy of the GCG and BM inside the 
horizon plus the entropy of the horizon do not decrease with time (Wang et al. 2006). 

The location of the apparent horizon f a in the FRW universe according to Cai et al. (2009) 
is obtained as 

f A = H- 1 (i + n k )- 1/2 . (14) 

For k = 0, the apparent horizon is same as the Hubble horizon, i.e. f\ = H~ l . 

Following Cai and Kim (2005), the associated Hawking temperature on the apparent horizon 
is given by 

^-tkk-mz)- (15) 

where 2 m A < ensure that the temperature is positive. 

Taking time derivative of both sides of (14) and using Eqs. (4), (5), (6), (7), (8), (11) and 
(12), one can rewrite Eq. (15) in real dimension form as 

TA = ^8vr(l+\)V2( 1 + ^- 3 ^Ch), 

= t~~ "j tz H n u/2 (l + n fc -g), (16) 
fc B 4vr(l + £l fc ) 1/2 

where h, c and &b are the reduced Planck constant, speed of light and Boltzmann constant, 
respectively. Now taking Hq = 72 Km s _1 Mpc _1 for the present time, one can estimate the 
Hawking temperature on the apparent horizon as 7a ~ ~ lO" 22 K. 

Following Maartens (1996), the BM (non-relativistic matter) temperature in the universe 
scales as T oc a -2 , so like Gong et al. (2007b) we assume here that the GCG temperature has 
a similar behavior T oc a~ n to avoid the negative entropy problem, where n is an arbitrary 
constant. It is not necessary to take n = 2 to ensure that the GCG is in equilibrium with the 
BM, since their dispersion relations could be completely different (see Lima and Alcaniz 2004; 
Gong et al. 2007a). Since the usual BM temperature in the universe decreases as the universe 
expands, we expect that the GCG temperature also preserves this property. Santos et al. (2006) 
showed that the GCG temperature in the adiabatic evolution from a dust-like to a de Sitter 
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cosmological model remains in the range < T < T*, where T* ~ 10 32 K, the temperature of 
the Planck era, is the maximum temperature of the GCG when it fills small volumes. Therefore 
the temperature of the universe filled with the GCG cools down as expected. This shows that 
the temperature behaviour of the GCG is different from that of the DE with constant EoS 
parameter u. The temperature of the GCG decreases instead of increasing. This tells us that 
the increasing or decreasing behaviour of the temperature of the universe dominated by DE 
is model dependent; it is not a general property associated with the DE (Gong et al. 2007a). 
A basic difficulty, however, is that the present-day DE temperature has not been measured. 
Lima and Alcaniz (2004) using a very naive estimate obtained the present value of the DE 
temperature as T^ E ~ 10 -6 K in a non-flat FRW universe and in the absence of interaction 
with BM. Zhou et al. (2009) showed that the DM temperature in the absence of interaction 
with the DE behaves same as the BM as Tdm ^ a <1 • They estimated T^ M ~ 10~ 7 K for the 
present time. Since the temperatures of the DE and the BM at the present time differs from 
that of the horizon T Ao ~ 10 -22 K, the systems must interact for some length of time before 
they can attain thermal equilibrium. Although in this case the local equilibrium hypothesis may 
no longer hold (Das et al. 2002, Wang et al. 2008; Pavon and Wang 2009; Zhou et al. 2009), 
Karami and Ghaffari (2010) showed that the contribution of the heat flow between the horizon 
and the fluid in the GSL in non-equilibrium thermodynamics is very small, O(10™ 7 ). Therefore 
the equilibrium thermodynamics is still preserved. 

The entropy of the universe containing the GCG and BM is given by Gibb's equation 
(Izquierdo and Pavon 2006a) 

TdS = dE + PdV, (17) 
where V = 47rf A /3 is the volume of the universe and T = T A . Also 

E=-^(j>ct + p b ), (18) 
3H 2 

P = P C h + -Pb = Pch = ^chPch = -3 — ^Chwch- (19) 

Taking time derivative of both sides of (17) and using Eqs. (4), (5), (6), (7), (8), (18) and (19), 
we obtain 

T A S = -4TrHr 3 A p Ch (l + u + w Ch )(l - , (20) 

where u = Pb/pch- 

Also the evolution of the geometric entropy on the apparent horizon 5a = Ttr\ (Izquierdo 
and Pavon 2006a) is obtained as 

T A 5 A = 47ri7f|p C h(l + u + w C h)(l " ^-) • (21) 
Finally, adding Eqs. (20) and (21) yields the GSL as 

TaSioi = 27rr A/ 9 C h(l + u + wch)^A, (22) 

where 5 to t = 5 + 5a is the total entropy. 

Taking time derivative of the apparent horizon (14) and using Eq. (4), one can obtain 

f A =^Hrlp Ch {l+u + ujch). (23) 
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Substituting Eq. (23) into Eq. (22) gives 

T A S tot = 8ir 2 Hr 5 A p 2 Ch (l + u + w C h) 2 > 0, (24) 

which shows that the GSL for the universe containing the interacting GCG with BM enclosed 
by the dynamical apparent horizon is always satisfied throughout the history of the universe for 
any spatial curvature, independently of the EoS parameter of the interacting GCG model. 
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